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CN ' Abstract 

We recently proposed a new approach to high energy nuclear scattering, which 
[ treats the initial stage of heavy ion collisions in a sophisticated way. We are able to 

■ calculate macroscopic quantities like energy density and velocity flow at the end of 

this initial stage, after the two nuclei having penetrated each other. In other words, we 
provide the initial conditions for a macroscopic treatment of the second stage of the 
collision. We address in particular the question of how to incorporate the soft compo- 
^ ll nent properly. We find almost perfect "Bjorken scaling": the rapidity coincides with the 

O , space-tune rapidity, whereas the transverse flow is practically zero. The distribution of 

i-^ ' the energy density in the transverse plane shows typically a very "bumpy" structure. 

X 1 Introduction 



Unfortunately there does not exist a single formalism able to account for a complete 
nucleus-nucleus collision. Rather we have to - at least for the moment - divide the 
reaction Into different stages (see fig. [l]) and try to understand them as well as possible. 
There is, first of all, the primary interaction when the two nuclei pass through each 
other. Since at very high energies the longitudinal size is, due to the gamma factor, 
almost zero (of the order 0.1 fm at RHIC), all the nucleons of the projectile interact 
with all the nucleons of the target instantaneously. In such a primary interaction many 
partons are created, which interact (in the pre-equilibrium stage) before reaching an 
equilibrium, referred to as quark-gluon plasma. The system then expands, passing via 
phase transition (or sudden crossover) into the hadron gas stage. The density decreases 
further till the collision rate is no longer large enough to maintain chemical equilibrium, 
but there are still hadronic interactions till finally the particles "freeze out", i.e. they 
continue their way without further interactions. 

The equilibrium stage is often treated macroscopically, by solving hydrodynamical 
equations. Here, usually very simplified initial conditions are used, like flat distribu- 
tions in space-time rapidity. This could be done much better by calculating the initial 
conditions for the hydrodynamical evolution on the basis of a realistic model for the 
primary interactions. 

We recently presented a completely new approach [|l], ^, |3[ for hadronic interactions 
and the initial stage of nuclear collisions, referred to as neXus, where we provide a 
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Figure 1: The different stages of heavy ion collisions. 



rigorous treatment of the multiple scattering aspect. Questions of energy conservation 
are clearly determined by the rules of field theory, both for cross section and particle 
production calculations, which Is not the case In all of the corresponding models used 
so far to calculate the Initial stage. In addition, we Introduced (currently only to leading 
order) so-called enhanced diagrams, responsible for screening and diffraction. It was not 
the Idea to create another model with some more features, but to provide a model which 
is theoretically consistent, and therefore much more realistic than all the approaches 
used before. We are therefore using neXus in order to determine macroscopic quantities 
after the first stage, when the two nuclei have traversed each other, such that these 
quantities may be used as initial conditions for a macroscopic treatment of the later 
stages of the collision. 

Calculating for example energy densities from a model like neXus or any other model 
for primary Interactions, is not trivial. We know the momenta of all the partons and 
we may calculate the energy density of the partonic system on a given hyper -surface 
(constant t), as shown in fig. ^. Here only the "resolvable" (or hard) partons are con- 



Figure 2: The naive method: only the resolvable partons are considered to determine the energy density on a 
hyper-surface (dots on hyperbola). 

sldered (dots on the hyperbola representing a hyper-surface). But this is certainly not 
the correct answer for the quantity of Interest, since there are many "unresolved" (or 
soft) partons around, which contribute to the energy density in a significant fashion. 
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The problem can be solved by using the string model, which is nothing but an attempt 
to treat the "soft partons" implicitly . The soft partons represent the string between 
"kinks", the latter ones representing the hard partons, see fig. ^. Thus we are able to 
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Figure 3: Between resolvable partons (thick lines) there are In addition unresolvable partons (thin lines), 
which are not treated explicitly. Their contribution is considered implicitly via strings. 

calculate energy densities in nucleus-nucleus collisions taking properly into account 
soft and hard (resolvable and unresolvable) contributions. 

The outline of this paper is as follows: we first review briefiy the basic elements of 
the neXus model, before we discuss in somewhat more detail the role of strings in this 
approach and the dynamics of strings. We then proceed to calculate energy densities 
for a given proper time, taking into account the hard and the soft partons. 

2 The NEXUS Model 

The most sophisticated approach to high energy hadronic interactions is the so-called 
Gribov-Regge theory Q. This is an effective field theory, which allows multiple inter- 
actions to happen "in parallel", with phenomenological objects called "Pomerons" rep- 
resenting elementary interactions . Using the general rules of field theory, one may 
express cross sections in terms of a couple of parameters characterizing the Pomeron. 
Interference terms are crucial, as they assure the unitarity of the theory. 

A big disadvantage is the fact that cross sections and particle production are not cal- 
culated consistently: the fact that energy needs to be shared between many Pomerons 
in case of multiple scattering is well taken into account when considering particle pro- 
duction (in particular in Monte Carlo applications), but not for cross sections f^. 

Another problem is the fact that at high energies, one also needs a consistent ap- 
proach to include both soft and hard processes. The latter ones are usually treated in 
the framework of the parton model, which only allows the calculation of inclusive cross 
sections. 

We recently presented a completely new approach [|l], ^, |3| for hadronic interactions 
and the initial stage of nuclear collisions, which is able to solve several of the above- 
mentioned problems. We provide a rigorous treatment of the multiple scattering aspect, 
such that questions of energy conservation are clearly determined by the rules of field 
theory, both for cross section and particle production calculations. In both (!) cases, 
energy is properly shared between the different interactions happening in parallel. This 
is the most important new aspect of our approach, which we consider a first necessary 
step to construct a consistent model for high energy nuclear scattering. 
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We first consider pp scattering. An elementary interaction is given as a sum of soft, 
semi-hard, and hard contributions: T2^2 = Tsoft + Tgomi + Thard. as discussed in detail 
in ref. |^]. We have a hard contribution Thaid. when the the first partons on both 
sides are valence quarks, a semi-hard contribution Tscmi. when at least on one side 
there is a sea quark (being emitted from a soft Pomeron), and finally we have a soft 
contribution, when there is no hard scattering at aU (see fig. T^ard is calculated 
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Figure 4: The soft elastic scattering amplitude Tgoft (left) , the hard elastic scattering amplitude Thard (middle) 
and one of the three contributions to the semi-hard elastic scattering amplitude T^cmi (right) . 



using the standard techniques of perturbative QCD, Tsoft is parameterized, and Tgomi is 
calculated as a convolution of Tsoft and Thard- We have a smooth transition from soft to 
hard physics: at low energies the soft contribution dominates, at high energies the hard 
and semi-hard ones, at intermediate energies (that is where experiments are performed 
presently) all contributions are important. 

Let us consider nucleus-nucleus (AB) scattering. The nucleus-nucleus scattering 
amplitude is defined by the sum of contributions of diagrams, corresponding to multi- 
ple elementary scattering processes between parton constituents of projectile and target 
nucleons. These elementary scatterings are the sum of soft, semi-hard, and hard contri- 
butions: T2^2 = Tsoft + Tsomi + Thard- A Corresponding relation holds for the inelastic am- 
plitude T2^x- We introduce "cut elementary diagrams" as being the sum over squared 
inelastic amplitudes, X)a'(^2^x)(T'2^x)*. which are graphically represented by vertical 
dashed lines, whereas the elastic amplitudes are represented by unbroken lines: 



This is very handy for treating the nuclear scattering model. We define the model via the 
elastic scattering amplitude Tab^ab which is assumed to consist of purely parallel ele- 
mentary interactions between partonic constituents, described by T2^2- The amplitude 
is therefore a sum of many terms. Having defined elastic scattering, inelastic scattering 
and particle production is practically given, if one employs a quantum mechanically self- 
consistent picture. Let us now consider inelastic scattering: one has of course the same 
parallel structure. Just some of the elementary interactions may be inelastic, some elas- 



tic. The inelastic amplitude being a sum over many terms - Tab^x = Yli '^ab^x ~ 

to be squared and summed over final states in order to get the inelastic cross section, 

which provides interference terms J2xi'^AB^x)i'^AB^x)* ■ These can be conveniently ex- 
pressed in terms of the cut and uncut elementary diagrams, as shown in fig. |^. One 
has to be careful about energy conservation: all the partonic constituents (lines) leav- 
ing a nucleon (blob) have to share the momentum of the nucleon. So, in the explicit 
formula one has an integration over momentum fractions of the partons, taking care of 
momentum conservation. This formula is the master formula of the approach, allowing 
calculations of cross sections as well as particle production. In the latter case, the mas- 
ter formula provides probability distributions for the momenta taken by the Pomerons 
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Figure 5: Example of a cut multiple scattering diagram, with cut (dashed Unes) and uncut (full lines) elemen- 
tary diagrams (Pomerons) . 

and the remnants. Avery detailed description with many applications and comparisons 
with data can be found in . 

So far we described only the basic version of the model. In reality we also consider 
triple Pomeron vertices to lowest order, as discussed in detail in ref. |^. We do not 
yet consider higher orders, nor do we consider the case where the two legs of the triple 
Pomeron are connected to different nuclei. All this is work in progress. 

3 Hadronic Structure of Cut Pomerons 

In order to develop our multiple scattering theory, we use a simple graphical represen- 
tation of a cut Pomeron, namely a thick vertical dashed line connecting the external 
legs representing nucleon components, as shown in fig. ^ This simple diagram hides 




Figure 6: Symbol representing a cut Pomeron. 

somewhat the fact that there is a complicated structure hidden in this Pomeron, and the 
purpose of this section is to discuss in particular the internal structure of the Pomeron. 

Let us start our discussion with the soft Pomeron. Based on Veneziano's topological 
expansion one may consider a soft Pomeron as a "cylinder", i.e. the sum of all possible 
QCD diagrams having a cylindrical topology, see fig. [?[ As discussed in detail in ||^] , the 
"nucleon components" mentioned earlier, representing the external legs of the diagram, 
are always quark-anti-quark pairs, indicated by a dashed line (anti-quark) and a full 
line (quark) in fig. |^. Important for the discussion of particle production are of course 
cut diagrams, therefore we show in fig. |7| a cut cylinder representing a cut Pomeron: the 
cut plane is shown as two vertical dotted lines. Let us consider the half-cylinder, for 
example, the one to the left of the cut, representing an inelastic amplitude. 

We may unfold this object in order to have a planar representation, as shown in fig. 
1^. Here, the dotted vertical lines indicate the cuts of the previous figure, and it is there 
where the hadronic final state hadrons appear. Lacking a theoretical understanding of 
this hadronic structure, we simply apply a phenomenological procedure, essentially a 
parameterization. We require the method to be as simple as possible, with a minimum of 
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Figure 7: Cut soft Pomeron represented as a cut cylinder. The grey areas represent unresolved partons. 
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Figure 8: Planar representation of a half-cylinder obtained from cutting a cylinder diagram (see fig. 



necessary parameters. A solution coming close to these demands is the so-called string 
model: each cut line is identified with a classical relativistic string. A Lorentz invariant 
string breaking procedure provides the transformation into a hadronic final state, see 
fig-i 

The phenomenological microscopic picture which stays behind this procedure was 
discussed in a number of reviews 0. H. the string end-point partons resulted from 
the interaction appear to be connected by a color field. With the partons flying apart, this 
color field is stretched into a tube, which finally breaks up giving rise to the production 
of hadrons and to the neutralization of the color field. 

We now consider a semi-hard Pomeron of the "sea-sea" type, where we have a hard 
pQCD process in the middle and a soft evolution at the end, see fig. 10. We generalize 
the picture introduced above for the soft Pomeron. Again, we assume a cylindrical 




Figure 9: The string model: each cut line (dotted vertical lines) represents a string, which decays into final 
state hadrons (circles). 
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Figure 10: A simple diagram contributing to the semi-hard Pomeron of the "sea-sea" type. 



structure. For the example of fig. 10, we have the picture shown in fig. 11: the shaded 
areas on the cylinder ends represent the soft Pomerons, whereas in the middle part we 
draw explicitly the gluon lines on the cylinder surface. We apply the same procedure 
as for the soft Pomeron: we cut the diagram and present a half-cylinder in a planar 



fashion, see fig. 11. We observe one difference compared to the soft case: there are three 





Figure 1 1 : Cylindrical representation of a contribution to the semi-hard Pomeron (left figure) and planar 
diagram representing the corresponding half-cylinder (right figure) . 

partons (dots) on each cut line: apart from the quark and the anti-quark at the end, we 
have a gluon in the middle. We again apply the string picture, but here we identify a 
cut line with a so-called kinky string, where the internal gluons correspond to internal 
kinks. The underlying microscopic picture will be presented by three color -connected 
partons - the gluon connected by the color field to the quark and to the anti-quark. The 
string model provides then a "parameterization" of hadron production, see fig. 12. The 
procedure described above can be easily generalized to the case of complicated parton 
ladders involving many gluons and quark-anti-quark pairs. One should note that the 
treatment of semi-hard Pomerons is Just a straightforward generalization of the string 
model for soft Pomerons, or one might see it the other way round: the soft string model 
is a natural limiting case of the kinky string procedure for semi-hard Pomerons. In a 
similar way one may treat Pomerons of valence type. 

The general picture should be clear from the above examples: in any case, no matter 
what type of Pomeron, unresolved soft partons play a very important role. In the string 
model, they represent the string pieces between the hard partons. In case of single 
Pomeron exchange in proton-proton scattering, particle production can be treated in a 
phenomenological fashion via the hadronization of two (in general kinky) strings. In nu- 
clear collisions, the situation is more complicated, since we have many Pomerons and 
consequently many strings - closely packed - which interact with each other. Never- 




Figure 12: The "kinky" string model: the cut line (vertical dotted line) corresponds to a kinky string, which 
decays Into hadrons (circles). 

theless, we can use the string picture to calculate energy densities, by considering the 
strings at an earlier stage, before they hadronize. 

It should be noted finally that particle production fi"om cut Pomerons is not the whole 
story. Cutting the complete diagram, one has as well to cut the projectile and target 
remnants, which may or not be excited. High mass excitations are as well considered 
as strings. 



4 Dynamics of Strings 

The string dynamics is derived fi-om the Nambu-Goto Lagrangian, which has been con- 
structed based on invariance arguments. The corresponding equation of motion for the 
string is a wave equation, with a solution [ll^. 111], fl^, 



r+t 



(1) 



for the four -vector X{r,t), having already assumed that the initial spatial extension of 
the string is zero. The quantity Xq represents the formation point of the string, which 
coincides with the position of the nucleon-nucleon interaction being at the origin of the 
string formation. The space-like variable r represents the position along the string for 
given time t, whereas the function g defines the initial velocity, 

g{r) = X{r,t)\t=o. (2) 

We will consider here a special class of strings, namely those with a piecewise constant 
function g, 

g{r)^Vk ior^^<r<—, l<k<n [3] 

K K 

for some integer n. We use for the string tension k = 1 GeV/fm. The set {E^.} is a 
partition of the interval [0, E], with E being the string energy, 

= Ea<Ei< ...<E,,_i<E„ = E, (4) 

and {ufc} represents n constant 4-vectors. Such strings are called kinky strings, with n 
being the number of kinks, and the n vectors Vk being called kink velocities. The function 
g must be symmetric and periodic, with the period 2E/n. This defines g everywhere, and 
eq. Q is the complete solution of the string equation, expressed in terms of the initial 
condition g. In the case of kinky strings the latter is expressed in terms of the kink 
velocities {vk} and the energy partition {Ek}. 
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What has all this to do with cut Pomerons? So far nothing, and to establish a link, 
we have to provide some mapping from the language of Pomerons and partons into the 
language of strings. We discussed earlier that a cut Pomeron may be identified with two 
sequences of partons of the type 



1-9-9 



9 - 9: 



(5) 



representing all the partons on a cut line. We identify such a sequence with a kinky 
string, by requiring 

parton = kink, (6) 

which means we identify the partons of the above sequence with the kinks of a kinky 
string, such that the partition of the energy is given by the parton energies. 



Ek = energy of parton k 
and the kink velocities are just the parton velocities, 

momentum of parton k 



Ek 



(7) 



(8) 



We consider massless partons, so that the energy is ec 



parton momentum. An example is shown on fig. [13 
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Figure 13: Partons of an e+e annihilation event with ^/s = 14 GeV In the z — x plane. The figure In the left 
upper corner shows the momenta In the — px plane. 

and an anti-quark, with 4 gluons in between - symbolically displayed in the first sub- 
figure, with a total cms energy of 14 GeV. One sees that the perturbative gluons play 
an important role in the beginning of the movement, and later from 2 GeV on, the 
longitudinal character dominates. A string breaks typically after 1 GcV/k, with n being 
the string tension, which gives much importance to the perturbative gluons. 
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5 Energy Densities from Strings 



Based on the string formalism, we are now able to calcvilate energy densities. As dis- 
cussed above, a string S* is a two-dimensional surface X{r, t) in Minkowski space, whose 
intersection with the surface 

t'-z'' = (9) 

defines uniquely a curve Sr, which we refer to as "String at given proper time r". Our 
aim is to calculate the energy density and the velocity flow at a given proper time r, in 
other words, along the curve Sr ■ 

Let us consider a small string piece A5t, in other words a segment of the string 
curve Sr between two neighboring points A and B, being sufficiently close so that the 
string piece may be considered as point-like. According to relativistic string theory, the 
four -momentum of this string piece in the lab system is given as 



Afl..= ; ,10, 



/;{ 

The partial derivatives can be expressed in terms of the Initial velocity g as 

dX{r,t) 1 



dt 2 
dX(r.f) _ 1 
Wr ^ 2 



[9{r + t)+g{r-t)\ (11) 
[g{r + t)-g{r-t)] . (12) 



This velocity function is known, in fact it is defined via the mapping of a system of n 
partons into string language: 

g{r) = Vk for < r < — , 1 < A; < n, (13) 

K K 

where Vk is the four -velocity of the k*^ parton divided by 7, and the difference Ek - Ek-i 
is its energy. So the four -momentum of a string segment can be easily expressed in 
terms of the original parton momenta. 
We use hj^erbolic coordinates, 

^ 1, 1 + z/t ^, 

qo=T = Vt^ - z\ qi=x, q2=y, = V = ^ log i_z/t ' ^ ' 

where z is considered to be the coordinate along the beam axis, and t is the time. The 

variables x and y are the transverse coordinates. For given r and 77, we define a frame 
via a Lorentz boost with boost rapidity rj. The four -momentum of the above-mentioned 
string piece is given as 

AP'' = A(;APi:b. (15) 

with the corresponding transformation tensor A. 

We are now going to calculate the energy momentum tensor T^" for a string piece in 
the frame F^. The general definition in kinetic theory is 



(16) 



where g is a position four -vector, and / the phase space density (particles per phase 
space volume) for a given time. For our point-like string piece we have in principle 

f{q,p)=6ip-AP)S{q-Q). (17) 
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Knowing that a string with zero width is a mathematical idealization, we introduce a 
Gaussian-type smearing function W{q), normalized as J W{q)4:TTq^dq = 1. So we define 



f{q,P) = S{p~ AP)W (jq- q||) 
with the norm in hyperbolic coordinates being given as 

\\q]\ = {q^f + iq2f + r^ {qzf . 



(18) 



(19) 



The energy momentum tensor for the complete string is then just the sum over all string 
pieces, which gives 



T'^'iq) 



E 



pO 



W 



q-Q 



string segments 

We define a local comoving frame via the four -velocity 

n^'{q) 



y/rU^iqyn^ 



with 



\q)^T°^iq). 



This allows then the calculation of the energy density in the comoving frame, 

eiq)^T'^^{q)u^{q)u,iq), 

and the flow velocity 

v{q) = u{q)/vP. 



(20) 

(21) 

(22) 

(23) 
(24) 



6 Results 

In the following, we show results for different high energy reactions. All the calculations 
are done at the proper time t = Ifm/c, and we suppress writing this variable in the 
following. We use always a width a ~ O.Sfm for smearing function W{q). We define the 
rapidity field yz{x, y, rf) as 



Vz{x,y,-q) = T^^og- -, 

2 l-Vz{x,y,ri) 



(25) 



where y, v^) is the z component of the velocity field. 

As a reference, we first show results for a single string without kinks, ha ving an 
energy of 50 GeV, as it may occur in electron-positron annihilation. In fig. 14, we 
plot the energy density e{x,y,r]) as a function of x and r/ for y = 0. As expected, the 
energy density is peaked around x = 0, whereas it is distributed evenly in 77 between 
limits defined by the energy of the string. Since there is a cylindrical symmetry with 
respect to the z axis, it is no surprise that the energy density as a function of x and y 
for r] = shows a narrow peak around a; = y = 0, as seen in fig. 15. The rapidity field 
yz{x, y, rj) as a function of x and r; for y ~ is shown in fig. 16, where we observe roughly 
yzix, y, T]) — rj. The other components of the velocity field are zero. 

Kinky strings show a very similar behavior. Due to the transverse momenta intro- 
duced via the kinks, the cylindrical symmetry is slightly distorted, and we observe as 
well very small but finite values for the transverse components of the velocity fields. Since 
the results are so close to the simple string discussed above, we do not show the figures 
here. 
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Figure 14: Energy density In the x — r] plane at y = for a single string of energy 50 GeV. 




Figure 15: Energy density in the x — y plane at ?? = for a single string of energy 50 GeV. 




Figure 16: Rapidity field In the x — rj plane at j/ = for a single string of energy 50 GeV. 
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Figure 17: Energy density In the x — ri plane at y = for a proton-proton collision at 50 GeV. 




Figure 18: Energy density In the x — y plane at ?? = for a proton-proton collision at 50 GeV. 




Figure 19: Rapidity field In the x — r) plane at y = for a proton-proton coUlslon at 50 GeV. 
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Next we consider proton-proton scattering at 50 GeV. Here we have usually several 
strings, the minimum number being two, plus two remnants (which may be strings in 
case of high mass excitation). In the example considered here, we have altogether four 
strings, with respective energies of 17.2 GeV, 9.3 GeV, 7.6 GeV and 15.8 GeV. Since the 
strings have different masses and are not sitting in the cms system, the energy density 
e{x, y, -q) as a function of x and -q for y = does not show such a flat behavior as in the 
case of a pure string, as seen in fig. 17, but we observe a peak at = 0, whose value 



is roughly 2.5 times bigger than the one for a string. However, the energy density as a 
function of x and y for 77 = shows as well a narrow peak around x = y = Q, as shown in 
fig. [1^, and the rapidity field yz{x, y, -q) is roughly equal to 77 , see fig. 19. 

Let us consider gold-gold collisions at RHIC (200 AGeV). In fig. |2C|, we show the 
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Figure 20: Energy density in the x — r) plane at ?/ = for a central gold-gold collision at 200 AGeV. 

energy density e{x,y,7]) as a function of x and ?/ for y = 0. For a fixed value of x, we 
observe a similar shape as for proton-proton scattering: a broad distribution with a 
smooth peak around zero. Of course, the magnitude is much bigger. Considering, how- 
ever, the variation with x for given ?/, we observed large fluctuations: pronounced peaks 
followed by deep valleys. If we regard the energy density as a function of x and y for 
77 = 0, as shown in fig. 21, we observe correspondingly several peaks overlaying the 




Figure 21: Energy density in the x — y plane at = for a central gold-gold collision at 200 AGeV. 

general roughly rotationally symmetric structure, which increases towards the center. 
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To investigate the origin of these fluctuations, we show in iig. 22 the distribution of the 
number of nucleons (projectile plus target) projected to the plane z = 0, in units nu- 




-8 -6 -4 -2 2 4 6 8 
x[fm] 



Figure 22: Projected nucleon density in the x — y plane for a central gold-gold coUlslon at 200 AGeV. 

cleons per fm^. From a spherically syrmnetric nuclear density, we expect a rotationally 
symmetric distribution of this projection, increasing towards the center {x — y — G). This 
is also what one observes, roughly. But looking more closely, we clearly observe large 
fluctuations with pronounced peaks. And even more, these peaks correspond exactly 
to the peaks in the energy density distribution, which proves that the fluctuations in 
energy density are due to geometrical fluctuations in the distribution of nucleons. We 




Figure 23: Rapidity field in the x — rj plane at j/ = for a central gold-gold collision at 200 AGeV. 



finally consider rapidity and velocity fields: in fig. 23 the rapidity field yz{x,y,if} and in 
fig. ^ the X component Vx{x, y, rj) of the velocity field, both as a function of x and rj for 
y = 0. We observe roughly yz = V and Vx close to zero. So the velocity field is practically 
purely longitudinal and shows the so-called "Bjorken scaling" (j/^ = r/). 



7 Summary 

We presented a new way to calculate macroscopic quantities like energy density or veloc- 
ity (or rapidity) fields at an early stage of a nucleus-nucleus collision at ultra-relativistic 
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Figure 24: Velocity field (x-component) in the x — r) plane at y = for a central gold-gold collision at 200 AGeV. 



energies. The calculation is based on a sophisticated treatment of the primary inter- 
actions, when the two nuclei are traversing each other, using the parton-based Gribov- 
Regge model neXus. The important point is an appropriate treatment of soft partons, 
which contribute substantially and which are usually completely neglected. In neXus, 
soft and hard physics are considered consistently: hard partons are treated explicitly 
based on pQCD, soft ones are included implicitly, using the string picture. This allows 
a quite reliable calculation of the above-mentioned macroscopic quantities. 

We analyzed single strings, proton-proton scattering and heavy ion collisions. In the 
latter case, we find almost perfect "Bjorken scaling": the rapidity coincides with the 
space-time rapidity rj, whereas the transverse flow is practically zero. This is often em- 
ployed as initial condition for hydrodynamical treatments. However, the -q dependence 
of the energy density does not show a well defined plateau corresponding to "boost in- 
variance". Furthermore, the distribution of the energy density in the transverse plane 
shows typically a very "bumpy" structure, which fluctuates considerably from event to 
event. 

This work has been funded in part by the IN2P3/CNRS (PICS 580) and the Russian 
Foundation of Basic Researches (RFBR-98-02-22024). H.J.D. acknowledges support 
from NASA grant number NAG-9246. 
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